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Complete expressions of time-domain gravitational waveforms for spinning binary inspirals via the 
post-Newtonian (PN) approximation would require determination of the phase, amplitude, inclina- 
tion angle, precession phase and spin vectors as well as the knowledge of the order coefficients for the 
PN expansion terms. These quantities are determined by solving simultaneously the spin-precession 
equations, the evolution equation for the Newtonian angular momentum, and the equation for the 
orbital frequency. For the spinning binaries with equal masses, determination of these quantities can 
be done fully analytically, by taking advantage of the equal mass symmetry, therefore by simplifying 
those equations to solve. We provide the analytical results through 1.5 PN order which includes 
spin-orbit interactions. 

Arun et. aZ (author?) [l| provides the "Ready-to-use" time-domain gravitational waveforms for spinning binary 
inspirals in Post-Newtonian expansion to 1.5 order. Their formulations, however, have yet to be further specified, 
depending on the configurations of binaries, such as the mass ratio and spin alignment with the Newtonian angular 
momentum. In order to make use of their theoretical waveforms in designing our wave templates via the IIR method, 
it is necessary to solve the equations governing these configurations and to fully specify the waveforms with all known 
parameters. 

The set of equations to solve are the following: 
(i) the spin-precession equations 



51 = J7i X Si, 

52 = f22 X S2, 



where at 1.5 PN order 



with 



M = Ml +M2, 
A'fiM2 

Ml - M2 



S = 



M 



(ii) the evolution equation for the Newtonian angular momentum 
where 

(iii) the equation for the orbital frequency 
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where 



Xs = ^(Xi+X2), (10) 

Xa = ^ iXl -X2) (11) 



with the normaUzed spin vectors 



X" = T?I' '^^l'^, (12) 



SO that \xn\ < 1 for objects that obey the Kerr bound on rotational angular momentum. 

I. SPECIFYING THE WAVEFORMS FOR EQUAL-MASS SPINNING BINARIES 

A. Configurations for tiie equal-mass case 

One interesting case might be a spinning binary with equal masses, where one can take advantage of the relatively 
simple spinning configurations due to mass symmetry, thus can analyze the spin effects on the binary such as precession 
and change in wave frequency more easily. For an equal-mass binary, one has 

M 

Ml - M2 = — , (13) 



therefore 
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4' 

5 = 0, (15) 
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and 



thus. 



-=1>2, (16) 



2S 



M2' 
2S 

X. = ^, (18) 



with the definitions 



S = S1+S2, (19) 
S = S1-S2. (20) 

With this simplification, one reduces the set of equations (i) , (ii) and (iii) above to the following: 
(i') the spin-precession equations 

51 = fix Si, (21) 

52 = nx S2, (22) 

with 

n^lM^'u^li^i^. (23) 



Further, we may combine Eqs. (|2ip and ([22|) via Eq. (|T9l) . 

S = r2 X S, (24) 
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which will be used throughout the rest of the analysis instead of Eqs. (|2ip and (|22p . 
(ii') the evolution equation for the Newtonian angular momentum 



via Eq. ([Tt 

(iii') the equation for the orbital frequency 
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B. Solving the configuration-equations in the Hmit S <^ L 



In order to solve the equations above effectively, one needs to prescribe the time-varying coordinates with respect 
to the fixed Cartesian coordinates, in which the Newtonian angular momentum becomes the same as that of a non- 
spinning binary. One such kind of the frame of coordinates may be written as 



(27) 



where l and a represent the angle of inclination (due to precession) and the phase of precession, respectively 
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cos t COS a cos l sin a 


— sini 






ey{t) 




— sin a cos a 













sin t cos a sin l sin a 


cos L 




ezo 



Ifj = e,it) e^o = Jo 




(f ) = cos t cos a e^o + cos i sin a — sini e^Q 
ey{t) = — sin cc + cos a By^ 
(t) = sint cos a c^q + sin i sin a e^o + cosi e^o 



The total angular momentum 

J = Jo = LAr + S (28) 

is assumed to be conserved (as radiation reaction is not included throughout our analysis) and is directed along the 
fixed z-axis in our analysis. We set Jq = e^Q and La? = ez(<), and via the relation of Eq. (P7| . Eq. (|28p may be 
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rewritten, 

Jo ^— sin /,ea;(i) + cos tl/AT^ = LnIjn + SS, (29) 
where Jq, £jv and S denote the magnitudes of Jq, L^v and S, respectively. From this one defines 



then finds easily 



C0S/3 = S-LAr = , (30) 



Ln + S COS (3 

COSL— . (31) 

Jo 



Also, from Eqs. (gHl) and ([3D]) we have 

= Ljf + 2LnS cos P + S^. (32) 
1. The angle of inclination 

According to Ref. (author?) fl'|, in the limit S ^ L the angle t can be considered a 0.5 PN correction, and taking 
advantage of this, we can reduce the lengthy expression for the GW polarizations to a much more compact form. 
Within this scheme, one finds out of Eqs. (|3ip and 



cost = 



yi + 2cos/3(^) + (^)' 
> 1 -Isin' 13 (^^y + 0[{S/Lr,f). (33) 



Now, in order to specify the angle i, we need to determine all the quantities involved in Eq. (j33p . First, the magnitude 
of the Newtonian angular momentum can be replaced by the leading order expression in ojorb, 

i^, = .A/^/ = 1m\-\ (34) 

where the latter expression is obtained using Eqs. (jS]) and ([T4|) . Next, from Eq. ([24]) . one finds that the magnitude of 
the total spin angular momentum S is constant (so is Xs) for our equal mass binary: 

i^(S'2) =S-S = S-(OxS) = 0. (35) 

From this fact and from Eq. ([50)1 together with Eqs. (|24p and (psj) . one also finds that the angle /3 is constant 

for the equal mass binary: 



d , ^ , A ^ S ■ Lat + S ■ Liv 4u S • S 
— (cos/3) = S-Ljv + S-LAr = ____ = o. (36) 



Then via Eqs. ([33]), ([34]), ([35]) and ([36]), one can finally specify 

cos . = 1 - ^^'^f'^ v' + O {v^) if S « L. (37) 

From this we may also infer 

45' sin /3 -f o ^ r fna^ 

t ^ j^^2 ^ if 5 < L. (38) 
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2. The orbital frequency 



The equation for the orbital frequency, Eq. p6)). can be integrated in a straightforward manner. First, we rewrite 
it using Eqs. (g]) and ((50)) . 



V = V 

5M 



Integrating this with respect to w, we obtain 
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where 



e = 



327r 376 5 cos ;3 



20M ' 



(39) 



(40) 



(41) 



and tc denotes the instance of coalescence, at which the frequency tends to infinity (the Post-Newtonian method 
breaks down well before this point) (author?) [2]. Now, one can invert Eq. (|40p and solve it for v, 



2 



i + i^e-V4 

4032 



TT 47 Scos(3 
10 ^ 120 M2 



ei/2 



(42) 



Via Eq. ([5]), one finds further 

'^"'■'^ - M = 8M® 



i + i^e-V4 

1344 



37r 47 5 cos ^ 
10 ^ 40 M2 



e-3/8 ^ 



ei/2 
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5. The precession frequency 



The precession frequency a also needs to be determined. Combining Eqs. p4p and (|25|) . and via Eqs. ([23]) . ([29 
and dMl) we find 



^-JoM^/'ujli^smLeyit) = -h f-Z^u^J^^' [ie,{t) + a sin Ley {t)] 



(44) 



Here, the ex{t) term being absent from the left-hand side can be justified by comparison of L and d on the right-hand 
side of Eq. (|44)) . To do so, we compare the ey{t) terms on the both sides first to find 
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Combining this with Eqs. ([5^ and ([M)) . we have 

7 , 



8M 



(45) 



(46) 



in the limit S <,L. Now, from Eqs. ([38]) and ([39]) we find 



32S'sin/3 
5Af3 ^ 



(47) 



Evidently, this quantity is much smaller than a therefore can be ignored in our analysis. By Eq. ([42]) we specify 
Eq. ^ further 
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4- The spin vectors 

As shown by Eq. ([35]) above, S, the magnitude of the total spin S is constant. So is Xs due to Eq. p7)) . 
First, the components of S in the basis {ex{t) , ey{t) , e2{t){= Lat)} can be found by solving Eq. ([24| . To do so, we 
insert 

S = 5^(*)e:,(t) + Sy'-''>ey{t) + S'^^'^Ln (49) 

into the equation. We have then, 

^ ■ 7nr2/^ 5/3' 

COS La ~ -gJVJ ' t^orb 

SvW ^ ^ (51) 

costd - ■gM'^/^LuJ^^ 

where the contribution from S'^(*)e^(t) + S'2'(*)ey(i) + S'^(*)LAr of S on the left-hand side of Eq. has been disregarded 

since its magnitude is much smaller than that of S'^(*)ei(t) + S^'^^^gyit) + S'^^^^Ln- From Eqs. §0^, ^ and (O 
above, we find 

S'^(*) = S- cos /3 = const. (52) 

Plugging Eqs. (gH), (gT]) and ^ into Eqs. ^ and and using Eq. (g]), we obtain 

S^W « ~S sin 13, (53) 
64 

SvW ^ —Ssm(3v^^0. (54) 
35 

Then via Eq. ^ [5'^(*), S'^'^*), fi-^^*)] transforms into [S"" , , 5^] in the basis {4o, eyo, e,o}- 

« S'(cos/? sin t — sin/? cos i) cos a « — S'sin/3 cosa, (55) 

S'^ « ^(cos/Jsint — sin/3cosi)sina « — S'sin/Jsina, (56) 

S^' « S'(sin/3sini + cos/3cosi) « S'cos/S. (57) 

Alternatively, we may express 

Xs = -XsSin/3cosa, (58) 

= — Xs sin/3sina, (59) 

Xf = Xs cos 13. (60) 

We can determine the components of S — Si — S2 in the basis {cxq, Cyo, Gzo} in a similar manner. First, one can 
show that the S, magnitude of S is also constant: so is Xa due to Eq. p^ . Subtracting Eq. from Eq. (PT|) . we 
have 

S^rjxS. (61) 

One finds trivially then, 

~{S^) =S-S^S-{nxS)^0. (62) 

In order to determine [S*''^*), 5^'^*), S'"^^*)] in the basis {e^it), ey(i), Lat}, similarly we insert 

S = 5"(*)e,(0 + Sy'^'^eyit) + S'^'^Ln (63) 

into Eq. (ICTI) and find 

■ 7»f2/3 5/3' ^ 
COSta — ' l^orb 

^^(*) = ^ (65) 

costd - ■gM'^/^ujJ^^ 
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Also, we define 

cos/3 = ^ (66) 

and have 

^^(*) = 5cos;3. (67) 

Now, substituting Eqs. (HZ]), dH]), (gT]) and (gTl) into Eqs. §^ and ((651), and using Eq. jS]), we find 

-^(t) _ gcos^sin^ ^ 
cos /3 

^ 64 5cos^sin/3^3^^^ 
35 cos p 



However, the magnitude S calculated by means of Eqs. (|67| . (j68)) and (|69|) shows 

(70) 

cos p 

and we easily find 

cos ^ cos /3 = const. (71) 

Then we may write 

S'^(t) ^ -5sin/3, (72) 
w 0, (73) 
« 5COS/3. (74) 

The rest of procedure to transform ^^^f*), 5^^*)] into [S'=,Sy, in the basis {e^o, Cyo, e^o} is the same as the 
above. We finally have 

5^ « -5 sin /3 cos a, (75) 

5^ « -5 sin /3 sin a, (76) 

S*" « Scos(3. (77) 



Or alternatively. 



Xa ~ -XaSin/3cosa, (78) 
Xa ~ -XaSin/Jsina, (79) 
Xl « XaCos/3. (80) 



C. Determination of the total phase and the amplitude factor 

1. The total phase 

Arun et. a/ (author?) [l| defines the orbital separation vector n(t), which is set to lie along ex{t) at initial time, 
i.e., n{t = 0) = ej;{t = 0), and rotates on the plane spanned by exit) and ey{t) by the cumulative angle $(t). Then 
one may write down the following two orthogonal vectors: 

h{t) ^ ex{t)cos<^{t) +ey{t) sin <^{t), (81) 
X{t) = -e;(i)sin$(t) +ej,(i)cos$(t). (82) 

We see that $(t) is the phase measuring how h{t) has rotated relative to the vector ex{t). For our precessing binary, 
however, ex{t) is itself rotating about Ljv(= ez(i)), which is associated with the angles l and a (see Eq. (|27I) V 
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Therefore the total rotation of h{t) about Ljv should be a combination of a rotation of n(t) with respect to the 
comoving basis etj(t) and ey{t), which is parametrized by ^(t) and a rotation of the basis due to the precession, which 
is parametrized by t and a. Ref.(author?) [l| shows this in the following way. In the basis {n, A, Ljv}, we introduce 
an "orbital-like" frequency Worb; which is defined as Woib = (v • A)/r. Then one may write 



■V — rh + rwoibA. 
Now, by means of Eqs. (I?7l) . ([5T|l and ((5^ one can show 

n = I $ + cos La) X ~ {i cos $ + sin i sin $d) L^r . 



(83) 



(84) 

However, by imposing Ljv = n x v/ |n x v| = n x n/ n x n , one finds that the term proportional to ILn in Eq. ((84|) 
must be zero. Thus, we have n = (v • A)A, and by identifying this with Eq. (|M)) via Eq. ([55]) and , we obtain 

i^oib = $ + cos id, (85) 

which may be now interpreted as the angular velocity with which n rotates about Ljv- The phase $(t) is then the 
integral 



$(<) = [ [Worb(i') - COS L{t')a{t')] dt'. 
Jq 



(86) 



By plugging Eqs. ([37]), (gll) and (j48]) into Eq. ^E^, and using Eqs. (gT]) and (|42|), one finally computes $(i) in the 
limit 5 < L: 



$(0 = -49^/^ 



^ U032 ^^96/ \ 4 ^^'^ 64 / ^ leVs 



(87) 



where pg = for the non-spinning case and pg = 1 for the spinning case, and S has been replaced by Xs via 

Eq. inD. 

Also, we need a separate expression for the phase of precession. Integrating Eq. with respect to t via Eq. (|¥T|l . 
and using Eq. ([T7)l, we obtain 



ait) = -Ps|e3/« 



(88) 



where ps = for the non-spinning case and ps = 1 for the spinning case. 



2. The amplitude factor 



Ref. (author?) [H gives the expressions for the waveform polarizations in the following form: 



rr(0) , fr(l/2) , rr(l/2.SO) „(1) 



„(3/2) „(3/2,SO) 



(l,SO) 



(89) 



In our equal- mass case (i/ = 1/4), the amplitude factor for each PN group will then be determined by means of 
Eq. JUl) 



i^(l/2 



i?(3/2 



M 



2Dl 8Dl 



e-i/4, 



M 



2Dl 16Dl 



M 



2Dl 32Dl 



e-V2 ("i + i^e-i/A 

V 1008 J ' 



M 



i + ^e-V4 



TT 47XsCOS^ I 3/8 

2 48 



(90) 
(91) 
(92) 
(93) 



2Dl MDl 4032 
where ps = for the non-spinning case and pg = 1 for the spinning case, and in the last equation S has been replaced 

2' 



by W'^Xs via Eq. 
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SUMMARY 



1. The total phase 



2435 
4032 



PS 



35 
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59xs cos (3 
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e-3/8 + O 



Ql/2 



2. The precession phase 



a(t) 



ei/4 



3. Inclination angle 



i « 2xs sin /3u if S" < i 



4. The spin vectors 



As 




— Xs sin ^ COS a 


Xs 




— Xs sin ^ sin a 






XsCOS^ 



Aa 




— Xa sin /3 COS a 


Xa 




— Xa sin /3 sin a 


Aa 




Xa cos^ 



5. The amplitude factors 



_p(l/2) 
^(3/2) 



2DI ~ MDl 



M 

M 

M 
32Dl 
M 



0-1/4 

e-3/8 
e-5/8 



487 
' 1008 
2435 



e-i/4 



4032 



e-i/4 



TT A7xs COS I3 \ 



Above, Ps — for the non-spinning case and pg — 1 for the spinning case. 
Xs = 2S/AP = constant and Xa = 2S/AP = constant; S = Si + S2, S = Si 

P = cos^i (t^) ^ (t^) ^ constant. 



S2 M = Ml + M2 
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